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HYCOM  2.2  (I) 


Maintain  all  features  of  HYCOM  2.1 

o  Orthogonal  curvilinear  grids 

o  Can  emulate  Z  or  Sigma  or  Sigma-Z  models 

o  It  is  “Arbitrary  Lagrangian-Eulerian”,  see: 
Adcroftand  Hallberg,  O.  Modelling  1 1  224-233. 

o  Explicit  support  for  1-D  and  2-D  domains 

o  KPP  or  Kraus-Turner  or 
Mellor-Yamada  2.5  or  Price-Weller-Pinkel 

o  Rivers  as  bogused  surface  precipitation 

o  Multiple  tracers 

o  Off-line  one-way  nesting 

o  Scalability  via  OpenMP  or  MPI  or  both 

o  Bit-for-bit  multi-cpu  reproducibility 

Special  halo  exchange  for  tripole  global  grid 

o  Arctic  dipole  patch  on  standard  Mercator  globe 
o  Logically  rectangular  domain 
o  Two  halves  of  top  edge  “fold”  together 
o  V-velocity  changes  sign  across  the  fold 


HYCOM  2.2  (Ha) 

Alternative  LeapFrog  barotropic  time  splitting 
o  Provided  by  SHOM 

o  Twice  as  expensive  as  standard  scheme 
o  Still  a  small  fraction  of  total  run  time 
o  Significantly  more  stable 
o  May  allow  2x  longer  baroclinic  time  step 

Alternative  scalar  advection  techniques 

o  Provided  by  Mohamed  Iskandarani 
o  Donor  Cell,  FCT  (2nd  and  4th  order),  MPDATA 
o  FCT2  replaces  MPDATA  as  standard  scheme 

Mixed  layer  changes 

o  GISS  mixed  layer  model 
o  Provided  by  Armando  Howard 
o  KPP  bottom  boundary  layer 
o  Provided  by  George  Halliwell 
o  KPP  tuning 


HYCOM  2.2  (lib) 

Initial  vertical  coordinate  changes 

o  Must  always  use  PCM  for  isopycnal  layers 

o  Vertical  remapping  used  PLM  for  fixed 
coordinate  layers 

o  Thin  deep  iso-pycnal  layers 

o  Stability  from  locally  referenced  potential  density 

o  Spatially  varying  layer  target  densities 
o  Different  isopycnal  layers  in  semi-enclosed  seas 

Recent  vertical  coordinate  changes 

o  hybrlx  only  active  below  “fixed  coordinate” 
surface  layers 

o  Major  re-write  of  HYBGEN  by 
George  Halliwell  and  Alan  Wallcraft 

o  Must  always  use  PCM  for  isopycnal  layers 
o  Vertical  remapping  uses  PLM  or  PPM  or 
WENO-like  PPM  (Laurent  Debreu)  for  fixed 
and  non-isopycnal  coordinate  layers 
o  More  layers  are  identified  as  non-isopycnal 
o  Updated  logic  for  two  layers 
(one  too  dense,  other  too  light) 


ANNUAL  MEAN  EQUATORIAL  PACIFIC 
GLBt0.72  HYCOM  VERSION  2.2.03 


u-velocity  zonal  sec.  O.OOn  mean:  4.004-  5.004  [06.0H] 


ANNUAL  MEAN  EQUATORIAL  PACIFIC 
GLBt0.72  HYCOM  VERSION  2.2.18 


u-velocity  zonal  sec.  O.OOn  mean:  4.004-  5.004  [22.3H] 


HYCOM  2.2  (lie) 

Equation  of  state  that  is  quadratic  in  salinity 

o  HYCOM  must  “invert”  the  equation  of  state 

o  tofsig(r,s)  and  sofsig(r,t) 

o  Traditional  version  is  cubic  in  T  and  linear  in  S 

o  Finding  the  root  of  a  cubic  is  expensive, 
but  exact 

o  Linear  in  S  is  not  accurate  at  low  salinity 
o  Optional  version  is  cubic  in  T  and  quadratic  in  S 
o  Coefficients  provided  by  Shan  Sun 
o  More  accurate  at  low  salinity 
•  Rivers,  Black  Sea,  Caspian  Sea 
o  Not  much  more  expensive 

Pade  equation  of  state 

o  Optional  Pade  version:  Pn/Qw 
o  P22  is  quadratic  in  T  and  S 
o  Qn  is  linear  in  T  and  S 
o  Developed  at  NCEP 
o  Only  SigmaO  at  present 
o  More  accurate  at  low  salinity 


HYCOM  2.2  (lid) 

Improved  thermobaricity 

o  Use  exact  formula  for  the  thermobaricity 

o  Eqn.  3  from  Sun  et  al.  1999 

o  No  single  reference  state  is  appropriate 
for  the  global  ocean 

o  Hallberg,  Ocean  Modelling,  8,  279-300 

o  Use  a  linear  combination  of  pressure  gradients 
from  two  out  of  three  reference  states 

o  Atlantic  (3°C,  35.0  psu) 

o  Arctic/Antarctic  (0°C,  34.5  psu) 
o  Mediterranean  (13°C,  38.5  psu) 
o  Most  locations  use  just  one  reference  state 

o  Linear  combinations  allow  smooth  transition 
between  states 

•  Do  this  in  shallow  water  if  possible 

o  In  deep  water,  constrain  the  T&S  used  for 
thermobaricity  to  be  close  to  the  reference  state 


HYCOM  2.2  (Ilia) 

Atmospheric  forcing  changes 

o  Option  to  input  ustar  fields 
o  Best  option  for  monthly  forcing 
o  Otherwise  calculated  from  wind  stress  or  speed 
o  Can  relax  to  observed  SST  fields 
o  Improved  COARE  3.0  bulk  exchange  coefficients 
o  Black-body  correction  to  longwave  flux 
o  Climatological  heat  flux  offset,  Qc 

Q  =  ( Qsw  —  Qiw )  +  (' Qi  +  Qs )  +  Qc 

o  Qc  is  constant  in  time 

•  Typically  based  on  the  model’s 
climatological  SST  error,  times 
(say)  -45  Wm-2/°C 

o  Analytic  diurnal  heat  flux  cycle 

o  Need  hourly  Qsw  for  good  cycle 

•  Have  3  or  6  hourly  (snapshots  or  averages) 
o  Input  daily-running  average  Qsw 

o  Apply  multiplicative  correction: 

•  clear-sky_now  /  clear-sky_daily-average 


SST  (degC) 


1  -D  HYCOM  at  OCEAN  STATION  PAPA 
NOGAPS  FORCING 


SST  in  2007  at  145W,50E 


Heat  Flux  (W/mA2)  SST  (degC) 


SST  in  2007  at  145W,50E 


days  since  Jan  1st,  2007 
Shortwave  Flux  in  2007  at  145W,50E 


days  since  Jan  1  st,  2007 


HYCOM  2.2  (Illb) 

Improved  support  for  rivers 

o  Still  bogused  surface  precipitation 
o  High  frequency  inter-annual  river  flow  allowed 
o  Add  it  to  atmospheric  precip,  off-line 
o  Instead  of  monthly  climatology, 
or  in-addition  to  it  (flow  anomalies) 

o  Better  control  of  low  salinity  profiles 
o  Option  for  mass  (vs  salinity)  flux 
o  Equation  of  state  that  is  quadratic  in  salinity 

Tidal  forcing 
o  Provided  by  NCEP 

o  Body  forcing  and  open  boundary  forcing 
o  Boundary  forcing  currently  for  “Flather”  ports 
o  v2.3:  Extend  it  to  Browning-Kreiss  nesting 
o  Body  forcing  for  8  largest  components 
o  SAL  treated  as  a  fraction  of  non-steric  SSH 
o  Tidal  drag  based  on  bottom  roughness 
o  Applied  to  near-bottom  tidal  velocity 
•  Use  a  lagged  25-hour  average  as  the 
non-tidal  velocity 


HYCOM  2.2  (IIIc) 

New  diagnostics  within  HYCOM 

o  Time-averaged  fields  (in  archive  files) 
o  Identical  to  off-line  mean  archives 
o  Daily  25-hour  average  for  tides 
o  No  in-line  capability  to  capture  variability 
o  Instantaneous  archives  still  available 
o  Sub-region  archive  files 
o  Example:  hourly  3-D  from  Global  1/12° 

•  3-4  small  regions  only 
o  One  file  per  involved  MPI  task  (entire  tile) 
o  Reconstruct  standard  regional  archive  files 
off-line 

o  Instantaneous  archives  still  available 
o  Synthetic  instrumentation 
o  Provided  by  George  Halliwell 
o  3-D  particle  tracking 
o  surface  and  constant  depth  drifters 
o  isopycnic  drifters 
o  fixed  instruments  and  moorings 


HYCOM  2.2  (Illd) 

I/O  optimizations 

o  Typically,  all  I/O  is  from  a  single  MPI  task 

o  I/O  can  be  a  bottleneck  when  running  on  many 
processors 

o  MPI-2  I/O  option 

•  Do  I/O  from  1st  MPI  task  in  each  row  of  tiles 
o  Faster  code  for  “endian”  conversion 

•  HYCOM  files  are  always  big-endian 

•  Intel  and  AMD  are  little-endian 

•  Do  the  conversion  in  parallel 
o  Sub-region  archive  files 

o  One  file  per  involved  MPI  task 
o  Can  be  much  faster  than  writing  a  full  archive 
o  Remove  density  from  restart  and  archive  files 
o  Less  I/O,  smaller  files 
o  Must  track  which  equation  of  state  is  used 


HYCOM  2.2  AND  SEA  ICE 


Finer  control  over  energy  loan  ice  model 

o  Melting  point  can  be  linear  in  salinity 
o  Set  ice  minimum  and  maximum  thickness 
o  Set  ice  vertical  temperature  gradient 
o  Or  get  ice  surface  temperature  from  Ta 
o  Made  compatible  with  coupled  sea-ice  approach 

Two-way  coupling  to  LANUs  CICE  sea  ice  model 

o  HYCOM  exports: 
o  SST,  SSS,  SSH 
o  Surface  Currents 
o  Available  Freeze/Melt  Heat  Flux 
o  CICE  exports: 
o  Ice  Concentration 
o  Ice-Ocean  Stress 

o  Actual  Freeze/Melt  Heat/Salt/Mass  Flux 
o  Solar  Radiation  at  Ice  Base 

o  Coupling  via  the  Earth  System  Modeling 
Framework 

o  Currently  for  non-global  domains  only 
o  v2.3:  ESMF  version  4  and  global  domains 


HYCOM  2.2  AND  CCSM 

Community  Climate  System  Model 
http  ://www.ccsm.ucar.edu/ 

o  Fully-coupled,  global  climate  model 
o  Sea-Ice:  CICE 
o  Ocean:  POP 

HYCOM  can  be  used  in  place  of  POP  in  CCSM3 

Uses  the  standard  HYCOM  source  code 

Subdirectory  CCSM3  used  to  hold  and  build  the 
CCSM3  version 

o  Some  source  code  files  are  specific  to  CCSM3 

o  HYCOM  ”.f”  files  are  renamed  ”.F”  to  simplify 
CCMS3  integration 

o  Macro  USE.CCSM3  for  CCSM3-specific  code 


HYCOM  2.2  (IV) 

Climatological  nesting  now  allowed 

o  Start  from  monthly  mean  outer  model  archive  files 
o  Allows  nested  runs  longer  than  the  outer  run 
o  But  with  less  accurate  boundary  state 
o  Probably  only  suitable  for  regional  nests 

Nesting  no  longer  requires  co-located  grids 

o  General  archive  to  archive  horizontal 
interpolation  (curvilinear) 

Hybrid  to  fixed  vertical  grid  remapper 

o  Allows  fixed-coordinate  nests  inside  hybrid 
coordinate  outer  domains 

o  HYCOM  to  (fixed-grid)  HYCOM 
o  HYCOM  to  NCOM 


HYCOM  2.2  (V) 

Enhanced  hycomproc  and  fieldproc 

o  NCAR-graphics  based 
o  Many  more  color  palette  options 
o  Can  read  in  an  arbitrary  palette 
o  Mark  locations,  and  draw  tracks,  on  plot 
o  Plot  diffusion  coefficients  and  tracers  (hycomproc) 
o  Overlay  vector  and  line-contours  (fieldproc) 
o  Vectors  can  optionally  follow  streamlines 

Added  fieldcell 

o  Like  fieldproc,  but  for  cell-array  (vs  contouring) 
o  Mark  locations  and  draw  tracks 
o  Overlay  line-contours 
o  Uses  NCAR’s  map  projections 

o  Typically  much  faster  than  fieldproc, 
but  can  leave  unfilled  cells 

o  Option  to  increase  resolution  of  input 
(bi-linear  interpolation) 


EXAMPLE  OF  FIELDPROC 
SSH,  SURFACE  CURRENTS,  AND  BATHYMETRY 
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FSD  (cm)  and  V.1  -  Jan  Year  20 


HYCOM  2.2  (VI) 

Diagnostic  fields  to  netCDF  and  other  file  formats 

o  Archive  fields  in  layer  space 
o  On  p-grid  (interpolated  velocity) 
o  3-D  archive  fields  interpolated  to  z-space 
o  On  p-grid,  or 

o  Sampled  at  stations  or  along  arbitrary  tracks 
o  3-D  archive  fields  sampled  on  iso-therms 
o  Meridional  stream-function  from  (mean)  3-D  archive 
o  In  logical  array  space  (rectilinear  grids) 
o  Binned  to  latitude  bands  (curvilinear  grids) 
o  Atmospheric  forcing  input  fields 
o  Time  axis  depends  on  “.b”  file  format 
o  Any  “.a”  file  with  the  right  “.b”  file  structure 
can  be  converted  to  netCDF 

o  Fields  binned  into  lon-lat  cells 


CANDIDATE  FEATURES  FOR  HYCOM  2.3 

Wind  drag  coefficient  based  on  model  SST 

Regional  tides  and  Browning-Kreiss  nesting 

Wetting  and  Drying 

Fully  region-independent 

o  Compile  once,  run  on  any  region  and  any 
number  of  processors 

o  Run-time  memory  allocation 
o  Might  reduce  performance  (fewer  compiler 
optimizations  available) 

o  Needed  for  full  ESMF  compliance 
Enhanced  support  for  ESMF 
o  HYCOM+CICE  on  global  tipole  grid 


